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On the Steady Motion of an Annular Mass of Rotating 

Liquid. 

By A. B. Basset, M.A., Cambridge, England. 



1. The recent investigations of Poincare* and Professor Q-. H. Darwinf 
have drawn attention to the problem of the figures of equilibrium of rotating 
masses of liquid ; and in the present paper it is proposed to consider the steady- 
motion of an annular mass of liquid whose cross-section is approximately circular, 
and which is rotating as a rigid body under the influence of its own attraction, 
about an axis through its centre of inertia which is perpendicular to the plane 
of its central line. 

This problem has to some extent been dealt with by Poincare, who has 
proved that such figures are possible forms of surfaces of equilibrium ; but the 
subject is capable of further development, and the object of this paper is to 
show how a solution may be obtained to any degree of approximation by the 
aid of the Toroidal Function analysis which has been so successfully employed 
by Mr. W. M. HicksJ in his investigations on circular vortex rings. 

2. In employing toroidal functions in problems such as the present, I have 
found it convenient to make use of a modified form of the methods introduced 
by Mr. Hicks. This method, together with many of the formulae required, will 
be explained in Chapter XII of my Treatise cm Hydrodynamics, but for the sake 
of completeness I shall proceed to give a preliminary sketch. 

Writing o = (a 2 + y z ) k , and starting with the dipolar transformation 

z + to = a tan — (£ + cyj) 



instead of the logarithmic form, and putting (7= cosh 77, c = cos£, k = e~ 



1 



* Acta Mathematica, Vol. VII, p. 259. t Phil. Trans., 1887, p. 379. 

X Phil. Trans., 1881, p. 609 ; 1884, p. 161 ; 1885, p. 725. 
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it has been shown by Mr. Hicks that the potential of an anchor ring or tore, 
which is composed of a mass of matter of constant density p, may be expressed 

in the forms y ={G + cf y °° ^ cos ^ 

*— « 'o 
at an external point, and* 

2 V~^ ™ 

where P n and Q n are the two zonal toroidal functions of degree n. 

From the formulae given by Mr. Hicksf it appears that P n and Q n respec- 
tively contain the factors 2h~ n + it and nh n+i , and we shall therefore find it con- 
venient to write 2PJc~ n+i and 7t^ m ^ +i for the functions which he denotes by 
P n and Q n . It will also be shown farther on that if (C + cf be expanded in a 
series of cosines of multiples of £, the coefficient of cos n% will be a rational and 
integral function of k multiplied by [2h)~ i . We shall therefore write 

V= (2J)*(C+ c y^A n P n (W* cos n£ (1) 

for the value of the potential at an external point, and 

V = - -|- Ttpr 2 + (26)*( ^ + c)* ]P 5„ Q n (&/&)• + * cos «g (2) 

for its value at an internal point, where b is the value of h at the surface of the 
tore. 

At the critical circle ri = oo , and therefore h=0. Now throughout the 
whole of the present paper the cross-section of the tore will be supposed to be 
small in comparison with its aperture, and as we shall only require to consider 
the values of the quantities in the neighbourhood of the tore, the following 
approximate values of P and Q and their differential coefficients with respect to 
Jc, which are denoted' by accents, will be sufficient for our purpose, viz : 

P Q = L+~^(L-1), P < = -h-i+±-k(L-\y 



p, 



= y (i + i-x*) f p> = b, p»=i, 

*p is supposed to be expressed in astronomical units. 
tPhil. Trans., 1884, equations (9) and (10). 



(3) 
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where L = log 4/ft ; also 



Q =l + -ft 2 , 



1 1 A 

Qo = ~2 ' Q = -g - , 

3 



*=!(! +-§-*•), «=!*• « ;,= ir 



, etc., 



5 
"16" 



ft=ir + 



etc. 



(4) 



3. We shall also require the expansion of r 2 , w 2 , and ((7+ e)*. 

Expansion of r 2 . 

t* __ C— c _ 1 — 2Jce 4- P 
If ~ C+e ~ 1 + 2ko + # 

= i _ 4fe. + sftV + 4ft 3 (c — 4c 3 ) — 16ft 4 (c 2 — 2c 4 ) , 

higher powers than ft 4 being neglected ; whence 

r 2 = a 2 {1 + 4ft 2 + 4ft 4 — 4 (ft + 2ft 3 ) cos £ + 4 (ft 2 + 2ft 4 ) cos 2£ 

— 47^ cos 3£ + 4ft 4 cos 4£ } . (5) 
Expansion of w 2 . 
jy _ (i_yy 
a 3 ~" (1 + 2ke + F} 2 

= 1 — 4ftc + 47? (3c 2 — 1) + 4ft 3 (5c — 8c 3 ) + ^ (1 — 9c 2 + 10c 4 ) , 
whence 
5 2 = a 2 1 1 + 2ft 2 + 2ft 4 — 4 (ft + ft 3 ) cos £ + 2 (3ft 2 + 2ft 4 ) cos 2£ 

— 8ft 3 cos 3£ 4- 10ft 4 cos 4£. (6) 
Expansion of (G + c) 1 - 

(2ft)*(C + c)* = (1 + 2&J + A?)* 

= l + &+-i^( 1 -c 2 )--i-^(c-c 3 )-^^ 4 (i-3c 2 +A c ^ 

whence 

(0+ c)*= (2ft)"* { i + -1ft 2 + -^-ft 4 + (ft- y ft 3 ) cos£ 

_ 2_#(i_ -L ftA cos 2£ + -1 ft 3 cos 3£ ^-ft 4 cos4£J. (7) 



4. Having obtained these preliminary results, we are in a position to con- 
sider the steady motion of the tore. 
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If £i be the angular velocity in steady motion, the condition to be satisfied 
at the surface of the tore is 

V + 4- & 0* 2 + f) = const. (8) 

We are not, however, at liberty to assume that the cross-section is an exact 
circle in steady motion, and we shall therefore suppose that its equation is 

e-' = h — I (1 + ft cos g + ft cos 2£ + ) (9) 

where b is the mean value of Tc at the surface, which is supposed to be small in 
comparison with a, the radius of the critical circle; and we shall also assume 
that ft, is a small quantity of the order b n . We must therefore first find the 
potential of an annular mass of matter of uniform density o whose cross-section 
is determined by (9). This is effected by assuming that the potentials at an 
external and internal point are respectively given by (1) and (2), and determin- 
ing the coefficients from the consideration that at the surface the values of 
J 7 " and V must differ by a constant, and that the values of dV/dk and dV'/dk 
must be equal. We shall thus obtain the values of the coefficients in terms of 
the fts. The resulting value of V and the surface value of g> 3 must then be sub- 
stituted in (8), and the coefficients of the cosines of multiples of £ equated to 
zero. This will determine the values of the fts. 

The preceding method will enable us to determine the values of the fts to any 
degree of approximation that may be desired, but in order to avoid unnecessary 
complication, the investigation will be confined to the determination of the first 
term of ft, which is sufficient to prove the existence of annular figures of equi- 
librium. From the course of the work it is evident that a higher approximation 
could be obtained with some additional labour. 

We shall find that B n is of the order b n , but that A n is of the order b n+ *, 
and we shall commence with the determination of V and dV'/dk. In calcu- 
lating the surface value of the former quantity, it will be unnecessary to proceed 
farther than the term involving cos 2£, or to include terms of a higher order 
than the second ; but in calculating d V'/dk it will be necessary, previously to 
performing the differentiation, to retain the term cos 3£ together with quantities 
of the third order in the coefficients, since the terms of the third order reduce 
upon differentiation to terms of the second order. 

5. Calculation of V. 

From (2), (5) and (7) we obtain 

V = — 4- »pa" { 1 + 4F— 4 (Je+2k 3 ) cos £ + 4*» cos 2£ — W cos 3£} + GH', (10) 
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where 
G = 1 + \# + (ft - —^ cosg - ~tf cos 2g + i-ft 3 cos 3£, (11) 

,5" = 5 g + 5 X & (ft/6) cos I + £,& (ft/6) 2 cos 2£ + 5 3 <2 3 (ft/6) 3 cos 3g. (12) 

Omitting terms of a higher order than the second, the surface values of the 
quantities are 

_ A Ttpr 2 = - -|- Ttpa 2 { 1 — 26ft + 46 3 — 46 cos g + (46 2 - 2&ft) cos 2£ f , 

(7=l+^6 2 + -^6 / 3 1 +6cosg + ^-(6/3 1 -i-6 2 )cos2?, (13) 

IT == 5, (l+ i-y) + -1^3!+ 4" A cob g + (^5^ + -|" B ») COB25 " 

Therefore the surface value of F' is 

F = const + (-|_ noa% + B 6 b + -|~ #i) cos £ 

+ {- 4 ^ 2 ( 4J - 2 W + T ^ x+ T 5l6 + T^ 2 1 cos 2 ^ (14) 

6. Calculation of d V /die . 

From (10), (11) and (12) we obtain 

— = 1 Ttpa 2 { 2ft — (1 + 67^ 2 ) cos g + 2ft cos 2g — 3ft 2 cos 3£ } 

NOW ^ _ 1 7 . , /-. 3 lA ™„ K 1 7, ™* 9.S -L- -1 

Also 



- _L jfe + (i _ 1^) cosg- Aft cos 2g + -g-7^ 2 cos 3g. (16) 



fl/ = Bo (l + A ft 2 ) + \ B, (l + A ft 2 )(ft/6) cos £ + A 5 2 (ft/6) 2 cos 2£ 

+ -A- J B 3 (ft/6) 3 cos3g. 



Therefore 



+ -A- 5 2 (ft/6) 2 } cosg + (~4^ + T ^'Z 6 ) C ° S 2 ^' (1?) 
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the term involving cos 3£ being omitted, as it is not required. Also 
-^ = \ BJc + \ BJr 1 (l + 4 #) cos I + 4" ^ & ~ % cos 2 £ 



1 ^ 
+ — B s b~W cos3£, 



therefore 
dW 



W 



+ 4 w^} cos s + (4- -w + 4" ^/ 63 ) c ° s ^ • ( is ) 

Substituting from (17) and (18) in (15), we obtain 
-jf — 3- "pa** + -Bo* + -J B i h l h 

+|4^ 2 (i+6^)+5 (i+4^)+4 5i ( 6 " 1+ 4 w ) + ^" jBa(W i co ^ 

+ J - -H *pa% -±BJc + ~ BJc/b + -| AW } co s 2 £ • 9 ) 

Putting for A; its value from (9), the surface value is 
^ = (_ ^ ^pa» 6 + J5 0& + -i ^)(1 + ft cos |) 

+ {4 7 tpa 3 (l + 66 2 ) + J5 (l + 4^)+^A&- 1 (l + 4 & + 4' jB2 } CO ^ 

+ {- -^ ^ - 4 ^ + 4 5i + 4" B * /h } (t & cos * + cos 2 0- (20) 

From (19) or (20) it appears that the most important terms ofdV'/dk are 
of zero order; and when we have calculated the value of dV/dk, it will be found 
that this circumstance requires that A n should be of the order b n + % . 

7. Calculation of V. 

Putting 

H = A P + A x (b/k) P 1 cos Z + A,(b/kYP 2 cos 2£, (21) 

we have V= OH. 
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The surface value of G is given by (13); also by (3) at the surface, 

P (k) = P (b) + -i fl-ft cos£ + (4"# -&) cos 2£, 
(6/fc) Pi cos f = cos £ — — ft — — ft cos'2£, 
(&/&) 2 P a cos2£ = -|cos2£. 
The surface value of i?is therefore 

+ | A (-i-fl-ft) -4-4A + ^ A} c os 2£. (22) 

Hence the surface value of Fis 

P = const + \A (Lb — ft) + A 1 \ cos 

+ |A(4^-/? 2 -i-^i)+4^( & -ft) + X^} COS ^' (23) 
8 . Calculation of dV/dk. 

We have dV _ Q dH R dO 

dh dk dk 

Now from (21) and (3), 

« =A j_*-. + A*( £ _4)}_^->{ 1 + i-*»(Z-i)|cOB| 

1 / 3 \ 
+ — ^ x 6 ( L —J cos £ + etc. 

This has to be multiplied by G whose term of lowest order is unity. Now the 
constant terms in dV'jdh are of the first order, and therefore the constant terms 
in dV /dk must also be of the same order, whence A /b must be of the first 
order, and therefore A must be of the second order of small quantities. Simi- 
larly it can be shown that all the other J.'s are of an order two degrees higher 
than their indices. We can therefore considerably simplify the calculation, since 
we do not require to retain quantities of a higher order than the second. We 
thus obtain 

G~ = - AJk - (AJ/tf + A ) cos I ; 
also from (16), 

H-jr- = A P cos £ = AJj cos £ , 
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W ieUCe ~dk = ~ A °l h + ( A ° L ~~ An ~ A h l^) cos I ' 

and the value of this at the surface is 

-£-=- All + (-A ft/6 + A L - A - A/b) cos g ; (24) 

the term involving cos 2£ being of the third order, is omitted. 

9. We can now obtain the value of V in terms of the fts. From (23) it 
appears that the coefficients of cos £ and cos 2£ in V are of the third and fourth 
orders respectively ; and since these quantities are equal to the coefficients of 
the corresponding terms in V, it follows from (14) that to the lowest order 

= A Ttpa'b + B b + ~B 1 , (25) 

= - ~ Ttpa 2 (4b"- 26ft) + A B x (b + ft) + A B t . (26) 

Equating the constant terms in dV/dk and dV'/d7c, we obtain 
-~7tpa 2 6 + B b+~B 1= - AJb, 
whence by (25) J, = A npaW. (27) 

Since the coefficient of cos 2£ in dV/dk is of the third order, it follows 
from (20) that 

— lfLnpaW — ±-BJ? + \ BJ> + ±B 2 = o. (28) 

Equating coefficients of cos £ in dV/dk and dV'/dlc and using (25) and (28), 

we obtain 

A (ft/6 + L — 1) - -4/6 = - 87tpa 2 6ft + 16?tpa 2 6 2 + -|- B,V + A ^J + A- £, 

= — 87ipa 2 6ft + 18?lpa 2 6 
by (28), whence 

A= A ?tpa 2 6 2 (16ft + 4Z6 — 316). (29) 

Again from (6), 

A fl 2 w 2 = A nV (1 — 46 cos | +~ ) . 

Substituting this together with the value of Ffrom (23) in (8), and equating the 
coefficient of cos t, to zero, we obtain 

ft= JL&(31 + 3n 2 /7tp6 2 — SL), (30) 

which determines ft. 
24 
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10. If we take the radius of the critical circle as the unit of length, it 
might be thought that if we gave to o any small numerical value, the value of 
£l would be arbitrary, subject only to the condition that the resulting value of 
/?! is a small numerical quantity of the same order as b. This, however, is not 
the case, for we have tacitly assumed that the pressure does not become nega- 
tive at any point in the interior of the ring. If the pressure did become 
negative at any internal point, this would indicate the existence of a hollow 
space within the ring, and the preceding investigation would be no longer appli- 
cable, for the potential at any point in the substance of a ring which contained 
a hollow would involve the P functions as well as the Q functions ; and we must 
therefore find the condition that the pressure should not become negative within 
the ring. 

Since V is the attraction potential, the pressure j> is determined by the 
equation 

Let P be the pressure along the critical circle where h = 0, then from (10), 
(11) and (12), 

P/p = - \ itpa 9 + B, + i-flV + C. 

We have already found the conditions that the coefficients of cos £ , etc., 
should vanish in p , hence to determine the condition that p should vanish at the 
free surface, all that we require is the constant term in V. By §6, this is 

-~7tpa 2 (1 + 4b- 26ft) + B + -l&(l? & + -i-i?i) + ~L ^(b^ + \b^ 

= _-|*pa , (l + 65») + .Bo 

by (25); whence 

= __ JL*pa»(i + Qtf) + B Q + |nV(l + 2V- 2b(3 1 ) + G 

and P/p = 47tpa?b* - Q?a? (6 a - &&) . 

We must therefore have 

47tp — fl a + &PJI > . (31) 

Let JL = Il 2 /47tp, then substituting the value of & from (30), the condition 
becomes # _ (i_ £ _ JJ| W + y > . (32) 
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11. As a numerical example, let b = .1 , then 

ft = .0124+ 10X, 

and if we put Jl = .01, the left-hand side of (32) becomes .9224, and is there- 
fore positive, and therefore 

ft =.1124, 

£l 2 /47tp = .01 
are solutions of the problem. 

Since (30) may be written in the form 



A = i(§ s ' + *-4«). 



it follows that if b is very small, ft cannot be of the same order unless % is very 
small. 

It is also necessary to point out that the preceding analysis proceeds on the 
assumption that ft is of the order Z> 2 , and it is quite possible that if a series of 
values were assigned to b and % which satisfy (30) and (32), these values might 
not satisfy the condition that ft should be of the order & a , or the inequality 
corresponding to (32) which would be obtained by carrying the approximation 
one degree higher. This difficulty would not be avoided by calculating the 
value of ft , since the same difficulty would exist with regard to ft . It there- 
fore appears that although toroidal function analysis throws some light on the 
solution of the problem, it fails to give a perfectly satisfactory approximate 
solution. The preceding results indicate that for small angular velocities an 
annular figure exists whose cross-section is approximately circular, but for large 
angular velocities the cross-section would probably become highly elliptical, and 
the ring would become flattened; and that if this quantity increased beyond a 
certain limit, the ring would break up. 

July 30, 1888. 



